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Abstract 

Einstein's field equations in FRW space-times are coupled to the 
BGK equation in order to derive the stress energy tensor including 
dissipative effects up to second order in the thermodynamical forces. 
The space-time is assumed to be matter-dominated, but in a low den- 
sity regime for which a second order (Burnett) coefficient becomes 
relevant. Cosmological implications of the solutions, as well as the 
physical meaning of transport coefficients in an isotropic homogeneous 
universe are discussed. 

1 Introduction 

The physics of low density systems is relevant in the description of the late 
universe. Since a hydrodynamical, Navier-Stokes, description of a low den- 
sity fluid breaks down for small Knudsen numbers, it is desirable to go be- 
yond this regime and analyze the fluid with the tools that statistical physics 
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provide as the particle collision frequency decreases. One simple tool use- 
ful for this purpose is the BGK equation, which simplifies the formalism 
of the full Boltzmann equation and provides a good grasp about how the 
non-equilibrium system evolves. In particular, dissipative effects, and gener- 
alized transport coefficients are derivable from the BGK formalism. These 
coefficients can be incorporated to the stress-energy tensor that is, in turn, 
coupled to Einstein's field equations. 

The relevance of cosmological transport coefficients has already been 
widely recognized |3J , nevertheless, until now, the second order Burnett equa- 
tions |2j have never been applied to isotropic-homogeneous space-times, 
although the issue has been addressed in the context of special relativity Q 
0. One possible motivation for the study of the Burnett-Einstein solutions 
is the well-known relation between wave scattering processes, density fluc- 
tuations and dissipative effects |J J/J. Dynamic structure factors, possibly 
relevant for describing CMB anisotropics or distortions, can be derived form 
the hydrodynamical Burnett approach, as it has already been done for colli- 
sionless plasmas fl5|-||10|. In this work we analyze the solutions to the BGK 
equation taking a FRW metric as power series expansions in the relaxation 
time r. To second order, the solution is used to obtain the constitutive equa- 
tions for the stress tensor thus providing explicit expressions for the transport 
coefficients. The explicit use of these results in specific situations will be de- 
ferred for future applications. Here we will be concerned with their physical 
and general features. Special emphasis will be placed in the significance of 
a bulk viscosity in homogeneous and isotropic systems. Unfortunate use of 
language has given rise to an identification of a geometrical aspect of the 
universe with a conventional dissipative effect. This gives rise to ambigui- 
ties that must be clarified. To accomplish this task, the paper is divided as 
follows: section 2 reviews the BGK equation [|I]]] in FRW space-times and 
outlines the basic strategy of its solution up to second order by means of a 
Chapman- Enskog type expansion (Burnett regime). Section three, is dedi- 
cated to the analysis of the first and second order distribution functions, the 
establishment of expressions for the transport coefficients and the construc- 
tion of the Einstein's field equation with a FRW flat metric. Final remarks 
about the implications of the Burnett regime in cosmology are included in 
section four. 
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2 BGK equation in FRW space-times 



We start this section writing the standard form of the BGK equation in its 
relativistic version for a single component system: 



£f_ 

Dt 



f ~ / (0) 
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where / = f(x fM ,v fl ) is the non-equilibrium distribution function, /(°) is the 
equilibrium function, r is the collision time (mean free time) and is the 
absolute derivative with respect to time t. The LHS of equation ([!]) can be 
written as: 

Df df „ , df , 



9f 4 , ®f a ri4 



a/3 



(2) 

reduces to 
(3) 



Dt dx» " ' dv» " v dx? ' " " a/3 
Assuming a flat, homogeneous and isotropic space-time, Eq 

Af _ 

Here we wish to make sure that the ensuing consequences remain clear. Eq. 
d^) has no explicit contributions arising from spatial gradients due to the 
assumption mentioned above. The second term in the R.H.S. containing the 
non- vanishing Christoffel symbols has an intrinsic geometrical origin so that 
all effects arising from it must be associated precisely to such structure. Now, 
denoting a(ct) the scale factor, c the speed of light, E the total mechanical 
energy of a single particle of rest mass m , 7 the usual relativistic factor and 
g^ u the metric tensor, and taking into account the expressions: 
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we can get a simpler form of Eq. (||), namely, 



D f E df df Ida ! i 22 33 

Dt=^m +m °dE { aTt ){VV +VV +VV) 



(7) 



or, since v a v 



-c, 

Df E df 
Dt m c 2 dt 



+ m, 



df .9.. E 



dE K 3^m 2 c 2 
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The factor 9 = is precisely a consequence of the geometrical features 



appearing in Eq. ([$]), and it is usually identified with the divergence of 
the four-velocity. Since the spatial gradients vanish in a FRW space-time, 
emphasis should me made on the fact that the non-vanishing terms of this 
divergence come from the Christoffel symbols introduced in the covariant 
derivative, so that 9 has a strictly geometrical origin. The solution of Eq ([!]), 
taking into account the expression (|8|) can be approximated up to second 
order in r using the proposal: 



/ = / (0) + rf»> + r^/W + 
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-2 f (2) 
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Substituting Eq. (|) into Eq. (H) and equating the coefficients of equal 
powers in constant r, up to second order, we get: 
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Expressions fllO-ll) constitute our first and second order corrections to the 
equilibrium distribution of the BGK equation in a flat FRW space-time. For 
non-relativistic particles in a FRW metric, Eqs. ( |I~0"l - |il"| ) reduce to: 



f (i) 

J NR 
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(12) 
(13) 
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The equilibrium function /'°) is, in this case, the Maxwell-Boltzmann distri- 
bution: 

r(0) , m 3 _™£ 771 3 _ E_ 

f () = n( — — ) 2 e 2kT — n( — r^) 2 e kT 14 
J K 2irkT J K 27rkT J v ; 

where n is the particle density, m is the mass of one particle of the simple 
system, k is Boltzmann's constant and T is the temperature. 
If particles are conserved, then the continuity equation holds: 

+ = ° (15) 



and the first term at the RHS of Eq. ([TJ) can be written as: 

-6f {0) (16) 



0/<o) df (o) dn 



dt dn dt 

In Eq. ( p!6|) 0, usually identified with V • u has, as emphasized before, a 
purely geometrical origin; it is, indeed the surviving term of the covariant 
derivative of a four velocity v a a . Straightforward calculations can now be 
carried, so that: 



m = (i + Uj (17) 



f (2) 



5mu 2 2 m 2 -u 4 2 
9A;T ~ 27 k 2 T 2 ' + dt 



/(°) (18) 



These corrections obviously vanish for homogeneous, isotropic, non-expanding 
space-times. The thermodynamical force inherent to the deviations from the 
perfect distribution is, to first order, the expansion scalar 6. To second order 
(Burnett regime), the generalized thermodynamical forces expected, 
non-linear forces appearing in the linear constitutive relations. Here, due to 
the assumption on space-time, they turn out to be simply 9 2 and 6%. The 
effect of these second order expressions on the field equations will be the 
subject of the next section. 

3 Stress-tensor and field equation 

The Einstein field equations with vanishing cosmological constant are given 
by: 

G£ = /cI? (19) 
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where is the Einstein tensor, k is the coupling constant and is the 
stress-energy tensor. The tensor Tjf is related to the non-equilibrium distri- 
bution function, up to second order, by: 



~m [ / (( V^ dV + mr [ / (1 dV + mr 2 [ f {2) v»v v dV (20) 



Eq. (|20|) is simply the ordinary stress-energy tensor written in terms of 
averages of the molecular velocities. The hydrodynamical velocity vanishes 
in the comoving frame. In this case, the first term of the R.H.S. of Eq. (p0|) 
may be written as: 



£1 47r 
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(21) 



where p is the mass-energy density, P = nkT is the pressure, h£ = 5£ + 
and f/ M is the hydrodynamic velocity given by: 
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The RHS of Eq. ( pij ) corresponds to the stress tensor of a non-dissipative 
fluid. If we wish to include dissipation up to second order in r for non- 
relativistic particles, we can use expressions flT7| - |IS| ) . Using this information, 
and assuming that no dissipation in the time axes occur, we can write for 
the first order correction of the stress energy tensor: 



(i) 

V 



' lrrW\ /(0) , d 
3 kT J 



9 = - Vc 9K (23) 



Following the rules of linear irreversible thermodynamics [13|, the relation 



(i) 



ship between the induced current and the force 9 which again we empha- 

(i) 

size arises from the inherent geometrical aspect of space-time, is T^= —r] c 9h^. 
T) c , the transport coefficient is usually referred to in the literature as bulk vis- 
cosity. This is misleading, since 9 does not contain spatial gradients. For 
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this reason, we shall refer to it as the cosmological viscosity r] c , here explic- 
itly given by: 

1 m u 2 



rj c = AirTm J ^1 



/<°> u 4 du 



(24) 



3 kT / 

The second order contributions to the stress tensor can now be identified as: 
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rpfj, 
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Thus, two Burnett transport coefficients arise, one for each generalized ther- 

dt 



modynamical force 9 2 , 9^, namely: 
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and the field equations (|T9[) , in the Burnett regime, in flat FRW space-times, 
take the rather interesting form 



(28) 



= ^pc 2 



(29) 



These equations, clearly dependent on the mean free time r, yield the dynam- 
ics of simple fluid consisting of non-relativistic particles in flat FRW metric 
in the Burnett regime. Future work will be dedicated to the study of their 
solutions. 

The order of magnitude of these coefficients is readily calculated. Substi- 
tuting /(°) one obtains that 



r] c m 2nkTr 



B l 



T 



and 



B 2 » -nkTr 2 
3 



(30) 
(31) 



The cosmological viscosity rj c will be relevant when r is of order one. 
Since for the Boltzmann dilute regime r ~ 10~ 5 s — 10~ 8 s, and for such values 
densities of the order 10 -5 , 10~ 19 g/cm 3 should prevail, corresponding roughly 
to the threshold of the matter domination era in the universe. Under these 
circumstances, the Burnett corrections which contain r 2 ~ 10~ 16 s 2 would 
require densities of the order of magnitude of 10~ ll g/cm 3 , the density of 
matter right after nucleosynthesis, an era where radiation is still relevant 
fL5|. Further study of these consequences will be reported later. 



4 Discussion of the results 



In this paper we have used the BGK model to describe the non-relativistic 
kinetic properties of a dilute gas in a FRW metric implying an isotropic and 
homogeneous universe. Two important results stem out of this analysis. In 
the first place the concept itself of transport coefficients is clarified. There 
existing no spatial gradients in the velocity, there will be no shear effects 
in the gas. Thus, ordinary viscosities, shear and volume, are not present. 
However, the contributions to the four dimensional covariant derivative of the 
velocity arising only from the geometric structure of space-time, embedded in 
the Christoffel symbols are present. They give rise to fluxes T£ arising from 
the dynamical character of the fluid. In the case here studied, the force turns 



out to be the expansion scalar 9. As we have shown elsewhere [|14||, there is a 



shear-like viscosity arising from geometry. Regretfully the nomenclature used 
in that paper is inappropriate, the coefficients then referred to are precisely 
rj c and its "shear-like" counterpart which should be called something else, 
not shear viscosity. 

The second important result is that as Eq. (|23|)-(p9|) point out, there was 
a Burnett-like hydrodynamic regime in the evolution of the universe which 
reflects itself in the field equations themselves, as well as in the transport 
properties of matter. Nevertheless, it did so at times when matter was still 
coupled to radiation, so that it may me useful to add to this analysis the 
behavior of radiation. This subject is, at present, one of study. 
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